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Abstract 


We have observed new channels for r decays with an r/ in the final state. 
We study 3-prong tau decays, using the r/ —> 77 and rj —► 37 t° decay modes 
and 1-prong decays with two 7 r°’s using the 77 —> 77 channel. The measured 
branching fractions are B(t~ —> -k~tt~tt + i]u T ) = (3.4/q ( ) ± 0.6) x 10 ~ 4 and 
B(t~ — > 7r _ 27r°7/n T ) = (1.4 ± 0.6 ± 0.3) x 10 -4 . We observe clear evidence for 
/1 —> i]7T7T substructure and measure B(t~ — > /i 7 r _ n r ) = (5.8tf 4 ±1.8) x 10 -4 . 
We have also searched for 7/(958) production and obtain 90% CL upper limits 
B{t~ —* 7r _ 7/V r ) < 7.4 x 10 -5 and t?(r” —> 7r~7r°7/u T ) < 8.0 x 10~ 5 . 

PACS numbers: 13.35.Dx, 14.60.Fg 
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Tau decays with an rj meson in the final state provide important information about 
various hadronic symmetries and allow for a study of the resonant structure of the weak 
hadronic current. These decays are rare and their detection became possible only recently 
with the high statistics CLEO experiment at Cornell Electron Storage Ring(CESR). Two 
such decays with small branching fractions already have been observed: B(r~ —> 'n~n 0 rji , T ) = 
(0.17 ± 0.02 ± 0.02)% 0 and B[r~ -> K~rjv T ) = (0.026 ± 0.005 ± 0.005)% @ 0. Both 
channels also have been seen by the ALEPH group 0. All other tau decays involving rj 
mesons were expected to be severely suppressed. The decay r —> 37 t rju T can proceed through 
the axial-vector current and its branching fraction was predicted to be 1.2 x 10” 6 ||. 

In this Letter, we present the first observation of the tau decay r —► 37 vr]v T using three 
final states: 7 r _ 7 r + 7 r _ ?/ 77 where rj is reconstructed from the rj —> 77 decay; 7T~7T + 7T~7j 3ir o 
where rj is reconstructed from its rj —> 3n° decay; and 7r _ 27r°?7 77 where rj is reconstructed 
from the rj —> 77 decay and the remaining photons form 2 7 r°’s. In addition, for the first 
time, we observe r~ —> /i 7 r _ z/ T using the /1 —1 7jn + n~ decay mode. We also search for decays 
with rj'{ 958) using the rj' —> ijTi + Ti~ decay mode with rj —> 77 . 

We use data obtained by the CLEO II detector 0 at the CESR operating at a center of 
mass energy corresponding to the peak of the T(4S') resonance ( E cm =10.6 GeV) and 60 MeV 
below this energy. The data correspond to an integrated luminosity of 4.68 / 6” 1 and contain 
about 4.27 million t + t~ pairs. CLEO II is a general purpose solenoidal spectrometer. In 
addition to good quality tracking, its special feature is a 7800 crystal CsI(Tl) electromagnetic 
calorimeter that provides photon detection with high efficiency and good energy and angular 
resolution, which is essential for rj and 7T° reconstruction. 

We select events using the 1 vs 3 and 1 vs 1 charged track topologies and tag one of the 
tau decays with a single charged track in the drift chamber which is required to be identified 
as an electron, muon or hadron. The electron candidate is required to have momentum, p, 
greater than 0.5 GeV and energy deposition in the calorimeter, E, such that 0.9 < E/p < 1.1. 
If specific ionization ( dE/dx ) information is available, we veto the event if it is more than 
two standard deviations below the expected value. Muon candidates must penetrate at least 
3 absorption lengths of material for track momenta less than 2.0 GeV, and more than 5 
absorption lengths for momenta above 2.0 GeV. 

A hadron tag is a track not identified as an electron or muon and with momentum pointing 
to the barrel part of the calorimeter, | cos 0 | < 0.81, where 6 is the polar angle defined with 
respect to the beam direction. The invariant mass, including all photon candidates in the 
tag hemisphere, is required to be less than 1.2 GeV. I 11 addition to single pions, this tag 
recovers unidentified electrons and muons and a large fraction of t~ —» p~u T decays. 

The second tau - representing the signal candidate - is reconstructed from its decays into 
7 T _ 7 r + 7 T _ rj 11 , 7 T _ 7 r + 7 T _ 7 / 37r o and 7r _ 27r°7 77 final states. We assume that all charged tracks are 
pions since there is very little phase space for decays in which one of the tracks is a kaon. 
The dE/dx information is consistent with this assumption. 

Photons are identified by isolated energy clusters in the calorimeter, separated from en¬ 
ergy deposits left by charged tracks and with photon-like lateral profiles of energy deposition. 
Photon candidates used for 7T° and rj reconstruction are required to be in the barrel part 
of the calorimeter and to satisfy |5/\ | < 10, where S* = (m iy — mx)/cr 77 (X = 7 r° or 
7j) and <t 77 is the 7 T° or 7 mass resolution(~12 MeV). Only photon pair combinations with 
—3.0 < S 77 < 2.0 are considered as signal candidates; those with larger values of IS/// are 
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used for sidebands. 

For lepton (hadron) tags, the lower energy photon used for r] reconstruction must have 
energy greater than 200 (250) MeV and the higher energy photon must have energy greater 
than 400 (700) MeV. Photons used to form 7 r°’s are required to have energies greater than 30 
MeV. In events for which more than one combination of photons passes all cuts, we choose 
the combination with the smallest y 2 for that signal hypothesis. 

We suppress e + e _ —> e + e _ ( 7 ) and e + e _ —> events by vetoing events with tracks 

which have energy greater than 85% of the beam energy. To remove background due to 2- 
plioton processes, we require the missing momentum vector of the event to be in the angular 
region |cosd| < 0.9. We suppress contributions from tau decays with a K$ in the final 
state by requiring that for all tracks, the impact parameter with respect to the interaction 
point must be less than 5 mm. Background from low multiplicity qq events and incompletely 
reconstructed tau events is minimized by rejecting events with additional isolated photons 
with an energy greater than 120 MeV. To further reduce qq background, we require the total 
invariant mass of the hadrons in the signal hemisphere to be less than the tau mass. 

For the sample we reduce qq and 2-photon backgrounds by requiring the 

event to have missing mass satisfying 0.1 < M miss /E cm < 0.5 and total transverse momen¬ 
tum greater than 0.3 GeV 0. Decays with AVs are additionally suppressed by requiring 
both 7 r + 7 r~ combinations to have a mass at least 15 MeV from the Ks mass. To suppress 
background from events with gamma conversions, we veto events with electron candidates 
in the signal hemisphere. 

We simulate tau signal and background events using the K0RALB generator and TAU0LA 
decay packages [§] (with some modifications discussed below) and measured tau branching 

m 


fractions 


GEANT 


is used for detector simulation. 


We find background associated with various QED processes to be negligible by using 
Monte Carlo (MC) simulations and independent data samples. To estimate qq background, 
we use independent data samples requiring the invariant mass of the tag hemisphere to be 
greater than 1.8 GeV. We select predominantly hadronic events satisfying the same topolog¬ 
ical and kinematic requirements on the signal hemisphere as described above, except for the 
tau mass cut. The normalization for this hadronic sample is obtained from a fit to the data 
in the region with signal invariant mass above 1.8 GeV. 

The distributions illustrating rj signal in all three analyses are shown in Fig. 1. Tau 
background contributes almost entirely through random 77 combinations while most of 
the coherent 7 background comes from qq events. Distributions of hadronic masses for 
events from the 7 / signal regions are shown in Fig. 2, where the events with masses above 
the tau mass are plotted as well. Signal regions for 7 —> 77 and rj —> 37 t° channels are 


3.0 < < 2.0 and |M(3vr° 




< 20 MeV respectively. 

We extract the number of r/’s by fitting the distributions shown in Fig. 1. In Table | we 
present results from the three analyses. Efficiencies shown in Tables I and II include tagging 
branching fractions. 

We estimate several sources of systematic errors. The major contributions are (for 
7r _ 7r + 7r _ ?7 77 , 7r _ 7T + 7r _ 73 ^ 0 , and n~2n°r]^ 1 samples): 7T° and 7 reconstruction efficiency (10%, 
10%, 4%); model dependence (5%, 10%, 10%); backgrounds (15%, 17%, 18%). The total 
systematic errors are 19%, 33%, and 20%, respectively. 

For the 7 r _ 7 r + 7 r _ ?/ 77 sample, there is enough data for a consistency check of B{t —> 
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FIG. 1. Distribution of rj mass for the (a) n~7r + n~rj 1 ^, (b) 7r~7T + Tr~rj 3n o and (c) 7r _ 27r°?7 77 
samples. The solid line is a fit to the data (squares). The tau and qq backgrounds are 
indicated by the dashed line and hatched area, respectively. Plots (a) and (c) are binned in 
units of S' 77 while (b) has a 10 MeV bin size. In (b), the energies and angles of each photon 
of the 7T° candidates have been constrained to 7 r° mass. 


TABLE I. The number of data signal events, efficiencies, and branching fractions for the three 
data analyses. The signal has been corrected for the background. 


Sample N)] ata e[%] £>(xl0 4 ) 

tt--tt+ <77--n 73 u+Dk FT q k+TTTn ; 


1 

+ 

1 

33 

3 

73.4l£“ 

6.3 

3.5l^±0.7 

1 

+ 

1 

-3 

CO 

o 

15-2 It! 

1.8 

3-il8;|±i.o 

1 

to 

o 

OS 

o 

o 

15-0±|;8 

2.5 

1.4lg;8±0.3 


7r _ 7r + 7r _ 77i/ r ) measurements among data samples selected with e, // and hadron tags. They 
are consistent with each other, with a y 2 of 5.1 for three degrees of freedom, corresponding 
to 16% confidence level. Combining 7r _ 7r + 7r“// 77 and -k~ti + ti~ r] 3n o results, we obtain B{t~ —» 
n~n + n~rju T ) = (3.41^5 ± 0.6) x 10 -4 . 

A 3nrj final state could proceed through a number of different resonances. In Fig. 3 
for events from the rj signal region we plot the Tirj vs ijim mass distributions, using the 
mass-constrained rj and 7r° momenta. The distributions show higher population density 
in the /i(1285)/ao(980) region, indicating the presence of the decay chain: r~ —> fi 
fi —> ao(980)7r, ao(980) —> rjn. For the 3-prong modes (Figs. 3(a) and (b)) there is an 
ambiguity in the choice of the same charge pions that results in four entries per event. In the 
case of r~ —> 7r~2n 0 rjh' T (Fig. 3(c)), there is only one fi —> rj^rr 0 combination and two rjn 0 
combinations. Since the kinematics of the fi —> ao7T° decay are such that the higher mass 
7jn 0 combination is the correct one about 90% of the time, we plot only the higher-mass ijn 0 
combination. 

We expect that more than 75% of all fi —> rjmr decays proceed through the a 0 (980)7r state 
['9j. To extract the number of the /i7r events we perform a binned maximum likelihood fit. 
We restrict the fit to the area shown in Fig. 3 to avoid the kinematically forbidden region and 
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FIG. 2. Hadronic mass spectra for the (a) 77 _ 77 + 77 _ ? 7 77 (after 77 sideband subtraction), (b) 
77 “77+77 “773^0, (c) 77 “ 277 °? 7 77 samples. The solid line is a fit to the data (squares). The tau 
and qq backgrounds are indicated by the dashed line and hatched area, respectively. 


we weight each event by the inverse of the number of entries. The fit function is the sum of 
a signal-MC ao vs fi distribution, a random 77 background shape obtained from non-signal 
tau MG and a constant background. For the signal MC we use the full decay chain with 
and a 0 (980) resonances 0 and include all 7/7777 and 7/77 mass combinations. The constant 
background accounts for a possible non-/! signal. To take into account the uncertainty in the 
random 77 background normalization, we perform a combined fit of the M( 7177 ) vs M( 77777 /) 
and S~ n data distributions with normalizations of the random 77 background constrained 
to be the same. From Monte Carlo studies we find that all three fits have confidence levels 
above 18%. We have found no sources of background which can contribute to the fi peak. 

I11 Table |TJ we show the fit results obtained for the different data samples. We have used 
£?(/i(1285) — > 777777) = 0.54 ± 0.15 |J and an isospin factor of 2/3 (1/3) for 7777+77“ ( 7777 ° 77 °). 
We include a systematic error of 28% to account for the uncertainty of the /1 —> 777777 decay 
rate. All other contributions to the total systematic error including different models of the 
fi —> 777777 decay are found to be much smaller. The total systematic error is 33% for all 
three channels. 

TABLE II. The numbers of signal events, efficiencies and branching fractions for the 
r“ —> fiTT~v r decays obtained from fits. 


Sample 

data 

£ [%] 

B{x 10“ 4 ) 

1 1 

+ + 

1 1 

CO 

3 

0 

36 . 3 / 9 / 

9-6±|f 

5.6 

1.4 

5.3/}/±1.8 

6 . 8 /|/± 2.2 

77 277°77 77 

8.4/1/ 

2.6 

6.6/1/T2.3 


The weighted average for all three channels is B(t~ —> f\T7~v T ') = (5.8//1 ± 1.8) x 
ICG 4 . Using results from the 77 “ 77 + 77 “ 7/ 77 channel we find the ratio B{r~ —» j\ i7~u T —> 
77 “ 7 r + 77 “ riu T ) /£>(r“ —> n~n + rjh> T ) = 0.55 ± 0.14. Here we take advantage of the fact that 
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FIG. 3. vs M(? 77 r 7 r) for the (a) 7 r _ 7 r + 7 r _ ? 7 77 (after 77 sideband subtraction), (b) 

TT~TT + TT~r] 3n o, (c) 7r _ 27r°?7 77 samples. Plot (d) is a sum of (a), (b) and (c) weighted as 0.25, 
0.25 and 1 respectively. All bin widths are 40 MeV. 


some systematic uncertainties are canceled in the ratio. It would appear that not all of the 
n~n + 7i~r] final state proceeds through an intermediate /j. We show in Fig. 4 the background- 
subtracted distribution of M n + n~ p) calculated for ti~tt + n~7]^ events with a 7r + n~7] mass 
of at least 36 MeV from the nominal _/j mass. This distribution as well as distributions of 
all other sub-mass projections are consistent with r~ —> ai(1260)“?p r —> tt~ p°rju T decay 
model. Since the observed excess has less than 2.0cr significance, we set an upper limit: 
B(t~ —■> ai(1260)~//z/ r —> 7T“ p°7)is T ) < 3.9 x 10~ 4 at 90% CL. 

The measured branching fraction for the decay r~ —> 7i~7r + T:~pi , T is more than two orders 
of magnitude larger than the value calculated by Pich || under the assumption that this 
decay is dominated by an rja\ intermediate state. In a recent calculation, Li |T2[ used similar 
assumptions to obtain B{r~ —> ai(1260)~7]u T —> ti~ p G r]v T ) = 2.93 x 10 -4 . In another paper, 
Li JTI| calculated B{t~ —> fin~u T ) = 2.9 x 10~ 4 , which is still somewhat smaller than the 
present measurement. 

The decay r —> 37has important implications for the phenomenology of the multi- 
pion r decays, especially for r —> 67T// t . Several authors [14-16] have used isospin relations 
17| to calculate the relative amounts of r~ —> 37r _ 27r + 7r°i/ r , r~ —> 7r _ 7r + 7r _ 37r°z7 T and t~ —> 
7 r _ 57r°z/ T , and claimed some discrepancies between the measured branching fractions and 
CVC predictions obtained from the e + e“ —> Gn measurements. 


It now appears that r 
contributions from the r — 


—> 3tt 27r + 7r ( V T and r —> 7 r 7 t + 7 t decays have large 

3'irr]i , T channel. Since this final state has opposite G parity to 
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FIG. 4. Distribution of ti~ 7i + 7i~ r t] invariant mass plotted for n~n + n ~events not associ¬ 
ated with t~ —> fin~u T decay. The Monte Carlo expectations for the decays t~ —> p° 7 r — 77 ^ 
and t~ —*• a!(1260) _ 7/z7 T —> 7T~7i + iT~r]i i T are shown as dashed and solid lines, respectively. 
Both Monte Carlo modes are generated via phase space with “V-A” factor Jll] . 


that of the direct 67 T decays and proceeds through an axial-vector current, its contribution 
must be subtracted before applying isospin relations or using CVC to compare with e + e“ 
annihilation data. 

We use the selected samples of r~ —> ti~ti + ti~ r]u T events to set upper limits on the 
t~ —> 7r - r/(958)ry and t~ —> 7r _ 7r°?7 / (958)z/ r decays. Each event must contain an rj candidate: 
0.51 GeV< M r] < 0.57 GeV. For the second decay, the remaining photons with E 1 > 30 
MeV are used for the 7T° reconstruction. I 11 the signal (sideband) region shown in 

Figs. 5(a,b) we find 2(1) and 0(1) events respectively. The event detection efficiencies are 
4.4% and 2.3%, with relative systematic uncertainties of 11% and 15.6% respectively. Using 
Poisson statistics and assuming a linear background distribution we obtain [[Dj: B{r~ —> 
n~rj'(958)u T ) < 7.4 x 10~ 5 and B{r~ -> 7r-7rV(958)i/ T ) < 8.0 x 10 ~ 5 . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun¬ 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftuug, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A.P. Sloan Foundation. 
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